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Abstract 

The Kolmogorov equation [|j with a forcing term is compared to experimen- 
tal measurements, in low temperature helium gas, in a range of microscale 
Reynolds numbers R\ between 120 and 1200. We show that the relation is 
accurately verified by the experiment (i.e. within ±3% relative error, over 
ranges of scales extending up to three decades). Two scales are extracted 
from the analysis, and revealed experimentally, one characterizing the ex- 
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ternal forcing, and the other, varying as i? A , defining the position of the 

maximum of the function — S^(r)/r, and for which a physical interpretation 

is offered. 
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Kolmogorov equation JT| is an exact relation between the longitudinal second order and 
third order structure functions, S2 (r) and S3 (r), valid for the ideal case of homogeneous 
isotropic turbulence. 5*2 (r) is linked to kinetic energy and S3 (r) to energy transfers, two cru- 
cial quantities characterizing fully developed turbulence. This relation is extensively used 
by the experimentalists to measure, from inertial range quantities, the mean dissipation 
rate e; no alternative method exists, in general, when the dissipative scales are unresolved, 
which typically happens at large Reynolds numbers. A restricted form of this equation, 
called "four-fifths law", is considered as one of the most important results in fully developed 
turbulence ||. The Kolmogorov equation was originally derived, after von Karman and 
Howarth [[| for freely decaying turbulence, and its adaptation to stationary forced turbu- 
lence, in a form suitable for a detailed comparison with experiment, was done by Novikov [[|; 
the corresponding equation, valid for scales well below an external forcing length Lf, reads: 



where r is the scale and v is the kinematic viscosity. Lf is an external scale, characterizing 
the forcing. We will call this equation "forced Kolmogorov equation" ; Lf is distinct from the 
integral length A, defined as the correlation length of the longitudinal velocity fluctuations. 
While the integral scale A is well documented, no systematic measurement of Lf has been 
reported yet. More generally, the extent to which forced Kolmogorov equation describes 
real turbulence is poorly known. The few investigations reported so far consider a truncated 
form of this equation — i.e. without the last term in Eq. ([!]) — , which, as will be shown 
in this Letter, is likely to be inaccurate for microscale Reynolds numbers R\ lower than 
approximately 1000. For larger R\, the truncated form is found compatible with the exper- 
iment, but sizeable deviations, on the order of 10 to 30%, are usually observed the 
existence of such deviations raises the issue as to whether the Kolmogorov equation should 
be amended to apply to real systems, and whether the fundamental concepts on which it 
relies — i.e. isotropic homogeneous turbulence — should be reassessed. Jeopardizing these 
issues, the results we present in this Letter show that the forced Kolmogorov equation de- 
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scribes the real world to a remarkable degree of accuracy, throughout the range of scales on 
which it is expected to apply. The analysis will further lead to single out two new scales for 
turbulence; one of them was introduced recently by Novikov 0, in a related context, but 
never observed. 

The set-up we use is the same as the one described in Refs. |]ID|-Q . The flow is confined 
in a cylinder, limited axially by disks equipped with blades, rotating in opposite directions, 
at approximately equal angular speeds. The working volume is a cylinder, 20 cm in diameter, 
and 13.1 cm in height. The cell is enclosed in a cylindrical vessel, in thermal contact with 
a liquid helium bath. The vessel is filled with helium gas, held at controlled pressure, 
and maintained between 4.2 and 6.5 K; the temperature is controlled with a long term 
stability better than 1 mK. Pressure and temperature are measured within 1% accuracy. 
The large scale structure of the flow is a confined circular mixing layer flQ |. Local velocity 



measurements are performed by using "hot" -wire anemometry. The sensors are made from 
a 7 yum thick carbon fiber, stretched across a rigid frame; a metallic layer covers the fiber 
everywhere except on a spot at the center, 7 /im long, which defines the active length of 



the probes. The time responses of the probes are analyzed, in some detail, in Ref . [O . We 
use here a probe located 4.7 cm from the mid plane of the system, and 6.5 cm from the 
cylinder axis; the speeds of the counter-rotating disks are finely tuned so as to maintain a 
local fluctuation rate close to 20%. We restrict our investigation to a range of R\ comprised 
between 120 and 1200. Here, R\ is defined by: 

* = (2) 

where u is the rms of the velocity fluctuations, and A is the Taylor microscale (based on the 
measurement of e discussed below). 

For all data files, more than 3 x 10 7 data points are recorded, ensuring comfortable 
convergence of the second and third moments. We check, for each file, that the velocity 
distribution is gaussian and we discard situations where the dissipative scale is not resolved. 
We also eliminate files for which the spectrum shows noise levels above 70 dB, or for which 



peaks (generally signalling the presence of a mechanical vibration) are visible in the inertial 
or dissipative ranges. This procedure leads to reject a little bit more than 50% of the files. 
For R\ > 1200, the noise level becomes prohibitive to ensure a reliable determination of 
S^r) and for R\ > 2300, we cease to resolve the dissipative scales. 

To investigate to what extent Eq. ([!]) agrees with the experiment, we write it in the 
following form: 

r + r dr 5* \ UL}) 1 } 

- following a procedure already used in Ref. |J — ; we will call the expression forming 
the left hand side J(r). All quantities on the left hand side can be accurately measured: 
the structure functions are obtained from the hot wire time series (we thus use Taylor's 
hypothesis to convert separation times into distances), and viscosity is known within ± 2% 
accuracy. The procedure thus consists in finding a best fit for the measured left hand side, 
by using the polynomial form given by the rhs of Eq. ([3|); in this calculation, two parameters 

- e and Lf — are free. The result is shown in Fig. [1] for R\ = 720. One finds the fit 
accurately reproduces the lhs of Eq. (|3]), within a range of scales covering two decades. The 
difference between the fit and experimental data is shown in the inset. The amplitudes of 
the deviations are below 3%, for r ranging between 8 and 900 r\. One could say that in 
this range, the forced Kolmogorov equation is verified within ± 3% accuracy. Outside this 
range, discrepancies are observed: below 8 rj, they are mainly due to noise which, although 
comfortably small for the usual measurements on turbulence, becomes here too large to be 
fully neglected. Above 900 rj, the discrepancies simply signal that the equation is no more 
applicable. 

Figure ^| collects results obtained for several R\, on Eq. (|I|), in the range 120-1200, i.e. 
one decade of variation in Rx- The points are measurements on the Kolmogorov function 
K{r) = —Ss/er and the full lines correspond to a determination of the rhs of Eq. ([I]), using 
best fit values for e and Lf, obtained by using the above procedure. As R\ increases, K{r) 
tends to form a plateau, in the inertial range, as expected from the Kolmogorov theory. 
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However, the trend, in terms of this parameter, is slow: inspecting the set of files we have 
for various R\ shows that, below R\ = 1000, there is no clear plateau, and one may ask 
to what extent an inertial range can be defined below this value. Anyhow, in all cases, 
the experiment confirms that the forced Kolmogorov relation accurately holds, the relative 
deviations between the theory and the experiment lying below 3%, for a set of records 
exploring three decades in scales, which is remarkable. 

We now turn to the analysis of the dependence on R\ of the characteristics of the curves 
of Fig. g. Those curves can be characterized by several quantities. One of them, the external 
scale Lf, is represented at various Reynolds numbers on Fig. [| There is a substantial scatter, 
but no systematic evolution with Rx is found, which indicates that Lf can be treated as 
a constant. We may thus consider that the effective forcing experienced by the flow is 
controlled by the flow geometry, which is physically acceptable. We estimate this scale as: 



We thus find Lf slightly smaller than the integral scale A (estimated to 2 cm in the present 
case), and one order of magnitude below the cell radius. The extent to which scales separate 
must be discussed by comparing Lf to rj. The fact that Lf is one order of magnitude smaller 
than the cell size implies that in practice one must achieve high Reynolds numbers to obtain 
scale separation and therefore fill conditions in which scaling behavior can be observed. The 
experiment also reveals an extended gap of scales, comprised between Lf and the cell size, 
for which turbulent fluctuations seem to escape from a theoretical description, based on 
simple assumptions. 

Another quantity of interest, useful for characterizing the plot of Fig. [| is the location of 
the maximum of K(r), which we call l s , for reasons which will appear later. By construction, 
this scale is well within the inertial range. l s is plotted against R\ in Fig. |3|; here again, 
there is a substantial scatter, but one finds a clear power law with Rx, in the form: 
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1.2 ±0.3 cm. 



(4) 



l s = (7.1±0.6)L f R x 



-0.57±0.04 



(5) 
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A last quantity of interest is the maximum value of K(r), whose evolution with R\ is 
displayed on Fig. |j. As expected, the maximum converges towards 4/5 as R\ increased. 
Nonetheless, the evolution is rather slow; the asymptotic regime is accurately reached (i.e. 
within 3%) only at R\ > 600. This observation, together with the preceding remarks on the 
formation of a clear plateau, underlines the fact that, in experimental systems, conditions for 
reaching the high Reynolds number limit, for the third order structure function, are difficult 
to achieve. The situation seems more favorable for the even order structure functions, and 
for the energy spectrum, the cross-over appearing at Lf being less pronounced, and hardly 
detectable if logarithmic scales are used. 

Simple characteristics of the maximum of the Kolmogorov function K(r) can be deduced 
from Eq. ([!]), by determining its approximate form, for scales well above rj, in a way similar 
to Ref. ||. Estimating the second order structure function S 2 by the expression: 

^(r) = co(er) 2//3 , where cq ~ 2, (6) 

and reinserting in Eq. (|3|), one gets: 

*<•■>= s-^r-Ms 1 < 7 > 

The maximum of K(r) is moreover found to be located at a scale: 

l s ~ LfR\~ 3 / 5 , (8) 
and the maximum of the function K(r) is calculated as: 
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with R\ ~ 30. The two results are well verified in the experiment: the value 3/5 we find is 
consistent with the experimental value 0.57 ± 0.04, and the expresssion for K max , plotted 
against -Ra in Fig. ^, agrees well with the experiment. We recall here that Eq. (|9|) applies 
only for moderate and large R\, i.e. well above 30, and should not be used for describing 
low Reynolds number turbulence, for which approximation (0) ceases to be valid. 



It is also interesting to reveal an equivalence between l s and another scale — which we 
temporary call I* — , introduced by Novikov || quite recently, and which was proposed to 
represent the size of vortex strings in fully developed turbulence. The equivalence between Z* 
and l s can be shown by reexpressing the quantity [called oi{r)] which controls, in the analysis 
of Ref. ||, the generation of vorticity correlations, and from which I* is defined. Assuming 
isotropy and homogeneity, one can derive, after some manipulations, the following exact 
relation between a(r) and S^r): 

1 ( d3 a d2 8 d 8\ PM 
a{r) = T4{ r ^ + 8 ^ + rTr--i) S ^- (10) 

In the analysis of Ref. @, I* is defined as the zero crossing of a(r); now, from Eq. (|IT)|). 
by using the approximation leading to Eq. ([/]), one can show that a(r) crosses zero axis at 
1.4 l s . Therefore, the two scales are equivalent, which justifies using a single notation for 
both. 



This remark further suggests a physical interpretation for l s . In a previous work [JT3 
internally coherent clusters of worms p^|15|| , whose size is proportional to ^ A ~ a70±010 have 
been found, a scaling close to l s . These clusters can be thought of as corresponding to the 
strings of Ref. since they define regions where velocity gradients (assumed to represent 



vorticity |14|]) are strongly correlated. We thus suggest here that l s provides a scale for the 
worm clusters in isotropic turbulence. 

To summarize, we have carried out a detailed, systematic comparison between a fun- 
damental relation of turbulence and the experiment. This study could be done by using 
low temperature helium, which allows to achieve highly controlled experimental conditions, 
while spanning an impressive range of Reynolds numbers. The analyses have been performed 
from single point measurements by using Taylor's hypothesis. We have shown that forced 
Kolmogorov equation applies to real flows, within a range of R\ lying between 120 and 
1200, with a remarkable degree of accuracy, estimated to ± 3% in relative magnitude, over 
three decades of scales. The experiment suggests that, as far as the third order structure 
function is concerned, the formation of an inertial range is slow. However, the results we 



found, whether close or far from the asymptotics, can be accurately interpreted by assum- 
ing an isotropic homogeneous turbulence state, which demonstrates the relevance of this 
approximation, used in almost all theoretical approaches to turbulence. The scales we infer 
from the analysis are observed for the first time, and we suggest here they may be useful to 
consider in order to characterize more completely experimental situations. 
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FIGURES 

FIG. 1. x: J(r) = —S^/r + 6uS' 2 /r compared to a best fit given by the rhs of Eq. (Q), for 
ii/\=720. In the inset, we show the difference between the best fit and the experiment; the full 
scale is ± 15%, and the dashed lines represent ± 3%. For this file, u/U = 20.9%, r] = 11.2 /zm, 
Lf = 964 r], l s = 147 r/ and e = 855 cm 2 s~ 3 . 

FIG. 2. The Kolmogorov function K{r) = —S^/er versus r/r], for different Reynolds number 
R\. The values of e are obtained by using best fits, as discussed in the text. v : -^A=120; O : -^A=300; 
A:i?^=1170. The solid lines show the expected curves, obtained from Eq. ([!]). 

FIG. 3. Scale Lf and the ratio l s /Lf, where is l s defined by the extremum of K(r), versus the 
Taylor Reynolds number R\. The solid line shows the best fit with an exponent -0.57 ± 0.04. 

FIG. 4. Evolution of the extremum of the Kolmogorov function, K(r) = —S^/er, versus the 
Taylor Reynolds number R\. The solid line shows the best fit given by Eq. (|9|), with R\ ~ 30. 
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